ABC as a flexible framework to estimate demography over space and time: some cons,many pros

The analysis of genetic variation to estimate demographic and historical parameters and to quantitatively compare alternative scenarios recently gained a powerful and flexible approach: the Approximate Bayesian Computation (ABC). The likelihood functions does not need to be theoretically specified, but posterior distributions can be approximated by simulation even assuming very complex population models including both natural and human‐induced processes. Prior information can be easily incorporated and the quality of the results can be analysed with rather limited additional effort. ABC is not a statistical analysis per se, but rather a statistical framework and any specific application is a sort of hybrid between a simulation and a data‐analysis study. Complete software packages performing the necessary steps under a set of models and for specific genetic markers are already available, but the flexibility of the method is better exploited combining different programs. Many questions relevant in ecology can be addressed using ABC, but adequate amount of time should be dedicated to decide among alternative options and to evaluate the results. In this paper we will describe and critically comment on the different steps of an ABC analysis, analyse some of the published applications of ABC and provide user guidelines.     

Size-dependent mortality in a Neotropical savanna tree: the role of height-related adjustments in hydraulic architecture and carbon allocation

Size-related changes in hydraulic architecture, carbon allocation and gas exchange of Sclerolobium paniculatum (Leguminosae), a dominant tree species in Neotropical savannas of central Brazil (Cerrado), were investigated to assess their potential role in the dieback of tall individuals. Trees greater than ∼6-m-tall exhibited more branch damage, larger numbers of dead individuals, higher wood density, greater leaf mass per area, lower leaf area to sapwood area ratio (LA/SA), lower stomatal conductance and lower net CO₂ assimilation than small trees. Stem-specific hydraulic conductivity decreased, while leaf-specific hydraulic conductivity remained nearly constant, with increasing tree size because of lower LA/SA in larger trees. Leaves were substantially more vulnerable to embolism than stems. Large trees had lower maximum leaf hydraulic conductance ( K _leaf) than small trees and all tree sizes exhibited lower K _leaf at midday than at dawn. These size-related adjustments in hydraulic architecture and carbon allocation apparently incurred a large physiological cost: large trees received a lower return in carbon gain from their investment in stem and leaf biomass compared with small trees. Additionally, large trees may experience more severe water deficits in dry years due to lower capacity for buffering the effects of hydraulic path-length and soil water deficits.     

Domestication of the Emblematic White Cheese-Making Fungus Penicillium camemberti and Its Diversification into Two Varieties

1. Download : Download high-res image (186KB)
2. Download : Download full-size image

     

Application of a photosystem II model for analysis of fluorescence induction curves in the 100 ns to 10 s time domain after excitation with a saturating light pulse

A mathematical model of photosystem II (PSII) events was used to analyze chlorophyll fluorescence transients in the time domain from 100 ns to 10 s after excitation with a saturating 10-ns flash, applied as a part of specialized illumination protocol, using preparations of a thermophilic strain of the unicellular green alga, Chlorella pyrenoidosa Chick (using both intact and diuron-treated cells). Analysis of simulation results has proven that particular attention should be given to flash-induced recombination processes, including nonradiative recombination in PSII, while subsequent charge transfer along the electron transport chain of thylakoid membrane can be adequately described by a single reaction of quinone reoxidation. The PSII model was extended by taking inhibition by diuron of the electron transport in the acceptor side of PSII into account, which allowed simulation of fluorescence induction curves observed in the presence of this inhibitor. The model parameters were determined (stromal pH, rate constants of nonradiative recombination, and the initial reduction state of the quinone pool) which provided adequate simulation of experimentally observed ratios of the maximal and initial fluorescence levels (F _m/F ₀).     

Consumer pressure, seed versus safe-site limitation, and plant population dynamics

Plants often suffer reductions in fecundity due to insect herbivory. Whether this loss of seeds has population-level consequences is much debated and often unknown. For many plants, particularly those with long-lived seedbanks, it is frequently asserted that herbivores have minimal impacts on plant abundance because safe-site availability rather than absolute seed number determines the magnitude of future plant recruitment and hence population abundance. However, empirical tests of this assertion are generally lacking and the interplay between herbivory, spatio-temporal variability in seed- or safe-site-limited recruitment, and seedbank dynamics is likely to be complex. Here we use a stochastic simulation model to explore how changes in the spatial and temporal frequency of seed-limited recruitment, the strength of density-dependent seedling survival, and longevity of seeds in the soil influence the population response to herbivory. Model output reveals several surprising results. First, given a seedbank, herbivores can have substantial effects on mean population abundance even if recruitment is primarily safe-site-limited in either time or space. Second, increasing seedbank longevity increases the population effects of herbivory, because annual reductions in seed input due to herbivory are accumulated in the seedbank. Third, population impacts of herbivory are robust even in the face of moderately strong density-dependent seedling mortality. These results imply that the conditions under which herbivores influence plant population dynamics may be more widespread than heretofore expected. Experiments are now needed to test these predictions. Received: 3 November 1999 / Accepted: 15 February 2000     

Adsorption and immobilisation of human insulin on graphene monoxide,silicon carbide and boron nitride nanosheets investigated by molecular dynamics simulation

The adsorption and immobilisation of human insulin onto the bio-compatible nanosheets including graphene monoxide, silicon carbide and boron nitride nanosheets were studied by molecular dynamics simulation at the temperature of 310 K. After equilibration, heating and 100 ns production molecular dynamic runs, it was found that the insulin was adsorbed and immobilised onto the considered surfaces in a native folded state. The structural parameters, including root-mean-square deviation and fluctuation, surface accessible solvent area, radius of gyration (R_g) and the distance between the centre of the mass of immobilised protein and the surface of the considered nanosheets, were measured, analysed and discussed. The energetics of the studied systems such as the interaction energy between protein and nanosheet was also measured and addressed. The discussions were centred on the structural and energetic parameters of the protein and nanosheets, including charge density, hydrophobicity, hydrophilicity and residue polarity. The results also showed that the active site of C-termini of chain B played an important role in the adsorption process and this could be helpful in the protection of insulin in its smart delivery and release applications.     

Structural variability of C3larvin toxin. Intrinsic dynamics of the α/β fold of the C3-like group of mono-ADP-ribosyltransferase toxins

C3larvin toxin is a new member of the C3 class of the mono-ADP-ribosyltransferase toxin family. The C3 toxins are known to covalently modify small G-proteins, e.g. RhoA, impairing their function, and serving as virulence factors for an offending pathogen. A full-length X-ray structure of C3larvin (2.3 Å) revealed that the characteristic mixed α/β fold consists of a central β-core flanked by two helical regions. Topologically, the protein can be separated into N and C lobes, each formed by a β-sheet and an α-motif, and connected by exposed loops involved in the recognition, binding, and catalysis of the toxin/enzyme, i.e. the ADP-ribosylation turn–turn and phosphate–nicotinamide PN loops. Herein, we provide two new C3larvin X-ray structures and present a systematic study of the toxin dynamics by first analyzing the experimental variability of the X-ray data-set followed by contrasting those results with theoretical predictions based on Elastic Network Models (GNM and ANM). We identify residues that participate in the stability of the N-lobe, putative hinges at loop residues, and energy-favored deformation vectors compatible with conformational changes of the key loops and 3D-subdomains (N/C-lobes), among the X-ray structures. We analyze a larger ensemble of known C3bot1 conformations and conclude that the characteristic ‘crab-claw’ movement may be driven by the main intrinsic modes of motion. Finally, via computational simulations, we identify harmonic and anharmonic fluctuations that might define the C3larvin ‘native state.’ Implications for docking protocols are derived.